Ultrafast Interference Imaging of Air in Splashing Dynamics 
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A drop impacting a solid surface with sufficient velocity will emit many small droplets creating 
a splash. However, splashing is completely suppressed if the surrounding gas pressure is lowered. 
The mechanism by which the gas affects splashing remains unknown. We use high-speed interfer- 
ence imaging to measure the air beneath all regions of a spreading viscous drop as well as optical 
absorption to measure the drop thickness. Although an initial air bubble is created on impact, no 
significant air layer persists until the time a splash is created. This suggests that splashing in our 
experimentally accessible range of viscosities is initiated at the edge of the drop as it encroaches 
into the surrounding gas. 

PACS numbers: 47.20Gv,47.20.Ma,47.55.D- 



When a liquid drop hits a surface, it may rebound pQ, 
spread smoothly, or shatter violently in a splash, as first 
photographed by Worthington [2 . Controlling whether a 
liquid splashes has important consequences in many ap- 
plications, including fuel dispersion in the automative in- 
dustry, splat formation in coating technologies, and pes- 
ticide application in agriculture [3j 0] . Liquid and surface 
properties obviously influence impact dynamics [5H9]; it 
is quite counter-intuitive however that lowering the ambi- 
ent pressure eliminates splashing altogether [T0HT3] . To 
measure transient air-layer dynamics, we develop a tech- 
nique that combines the high spatial precision of inter- 
ferometry (nm scale) with high time resolution (15/is). 

At impact, a small amount of air is trapped beneath 
the falling drop, creating a bubble [T4HT8] . Recent theo- 
retical work has suggested that this air pocket is linked 
to splashing dynamics [T9H2T] . In a sufficiently viscous 
liquid, splashing occurs at late times, several tenths of 
millisecond after impact [TTJ [12]. This temporal sep- 
aration between impact and splashing creates an ideal 
system to test whether the initial air pocket influences 
the later- time splashing dynamics. Using our interfer- 
ence technique, we find the initial air cavity dynamics to 
be consistent with theoretical predictions [20] , However, 
we find no significant air layer that persists beneath a 
spreading drop until the time of thin-sheet ejection — 
a necessary precursor to splashing in high-viscosity liq- 
uids [12 j. Thus, an underlying air layer is not responsible 
for splashing in this high- viscosity regime. 

The drops used in this study were mixtures of wa- 
ter and glycerol, with kinematic viscosities, between 
9 and 58 cSt. Over this range there is minimal varia- 
tion in surface tension, 7 (65-67 dyn/cm) and density, p 
(1.1-1.2 g/cm 3 ). Drops of uniform radius, ro=2.05 mm, 
were created using a syringe pump. For each impact, we 
used a fresh glass substrate (Fisherbrand coverslip) to 
prevent surface contamination. The impact velocity, ^0, 
was varied between 1.5 and 4.1 m/s by releasing drops 
from various heights within an acrylic tube that could 
be evacuated to varying ambient pressures, P, between 



2 and 102 kPa. Using a high-speed camera (Phantom 
vl2, Vision Research), we imaged drop impacts as shown 
in Fig. la. To determine the thickness of the liquid as 
a function of position and time, we measured the local 
optical absorption of a spreading drop of colored liquid 
(Brilliant Blue G dye in a v=9 cSst glycerol/ water solu- 
tion) as shown in Fig. lb; we converted the transmitted 
light intensity to liquid thickness by calibrating with a 
liquid wedge of known proportions. By modifying our 
setup as shown in Fig. lc, we also measured the thick- 
ness of any air layer underneath the spreading liquid us- 
ing interferometric high-speed imaging at speeds up to 
67,000 frames/second. We used a monochromatic LED 
(A=660 nm) with a small coherence length, ~ 10/im, as 
a light source, so that there would be no interference 
between the two sides of the glass substrate. Adding a 
small amount of dye to our liquid greatly minimized the 
reflected light from the upper liquid surface and elimi- 
nated any interference generated within the liquid itself. 

Fig. la shows a ^=15 cSt drop at atmospheric pressure: 
after spreading smoothly as a thick lamella for ~ 0.4ms, 
it ejects a thin liquid sheet that subsequently disinte- 
grates into smaller droplets — the splash. Using optical 
absorption, Fig. 16, we find the lamella edge to have a 
thickness of 106 ± Apm, while the ejected sheet is ten 
times thinner, only 10 db 2jim thick. This jump in thick- 
ness occurs over a lateral extent of only ~ 300/im. 

There are many distinct splashing regimes that dis- 
play different scalings and even qualitatively different 
behavior [3] [22]. For example, below ~ 3 cSt, sheet 
(i.e., corona) formation occurs within a few jis of im- 
pact, while above 3 cSt sheet ejection is delayed [TT] . 
While it is not at all clear that the instability is the same 
across all splashing regimes, it is nevertheless established 
that lowering the air pressure eliminates splashing in all 
cases [lOj [TT] . The higher viscosity liquids used in this 
study allow for a large separation in time and space be- 
tween the initial air layer entrapment and the creation 
of a splash. This allows us to test directly whether the 
initially trapped air layer persists to longer times to in- 
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FIG. 1: Images of a spreading drop after impact onto dry glass, (a) Side view showing thin-sheet ejection and splashing for 
a drop with v—Vo cSt, and wo=3.5m/s. (b) Bottom view and intensity profile of a dyed z/=9 cSt,i£o=3.3m/s drop used to 
meaure liquid thickness: the ejected sheet is approximately 1/10 the thickness of the lamella. Arrows indicate measured regions, 
(c) Schematic for reflected- light interference: a half-silvered mirror (Ml) directs light from a monochromatic LED onto the 
impacting drop from below; a standard mirror (M2) redirects reflected light into a high-speed camera, (d) Example time-series 
interference image of an impacting drop with v—2\ cSt and uo=3.5m/s. The lamella is uniform in intensity, indicating an 
optically flat region. The thin sheet is ejected in the third panel, (e) Magnified images of initial air cavity. This initial air 
cavity closes completely well before the thin sheet is ejected, (f) Radius of curvature of cavity normalized by drop radius, 
Tc/tq vs. time showing the cavity is slender (r c ~ ro), and gently flattens in time. Complete cavity closure occurs ^150 fis 
after impact, well before splashing occurs. 



fluence the splashing dynamics. 

Our interference technique determines the air-layer 
thickness beneath the drop as it spreads. Fig. Id shows 
an example interference image at several different times 
after impact: (i) just after impact, there is a small air 
cavity (panel 1) (ii) surrounding this tiny region, the 
spreading lamella is uniformly black indicating an op- 
tically flat surface (panel 2), and (Hi) underneath the 
ejected thin sheet (panels 3-5), interference fringes are 
widely spaced indicating a very shallow slope. 

The small cavity of air trapped under the impacting 
drop (Fig. le) has been shown to be present under vary- 
ing conditions, [T4HT8] including above and below the 



splashing threshold [12]. By measuring the interference 
fringes that are clearly observed inside the cavity, we can 
directly measure the cavity curvature as a function of 
time. The air cavity is quite fiat — the radius of cur- 
vature at the top of the cavity, r c , is comparable to ro- 
At impact, the overpressure at the edge of the cavity is 
predicted to be higher than at its center [20]. This sug- 
gests r c /ro should increase in time as the overpressure 
causes the cavity to flatten; this is consistent with the 
data shown in Fig. If. 

The air cavity persists even as the ambient pressure 
is decreased below the threshold value for sheet ejection, 
P S h [12] • However, the shape of the cavity strongly de- 
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FIG. 2: The radius of curvature of the entrapped air cavity 
normalized to drop radius, r c /ro, vs. ambient pressure P. In- 
set schematic defines cavity height, ft c , radius of curvature, r c , 
and lateral cavity radius, A. As P decreases, r c /ro increases 
and the cavity becomes flatter. The line shows the best fit to 
r c /r oc P _1 . 



pends on pressure. As shown in Fig. 2, the cavity flattens 
dramatically with decreasing pressure: r c /r oc P- 1 . To 
compare our curvature measurements with the theoreti- 
cal prediction [20] that cavity height, ft c oc P, we approx- 
imate the cavity as a thin spherical cap: 
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where 2A is the lateral extent of the cavity (see Fig. 2 
inset). Our measurement that r c /r Q oc P _1 thus corrob- 
orates the prediction that ft c oc P. We further note there 
is no transition in r c /ro at P sh , emphasizing that bubble 
entrapment does not appear to be related to the insta- 
bility producing the thin sheet. Moreover, we note that 
the cavity closes into a bubble well before sheet ejection 
occurs. We therefore conclude that the cavity dynamics 
is isolated in time and space from the edge of the spread- 
ing lamella where thin-sheet ejection and splashing take 
place. 

Underneath the lamella, outside the bright spot cre- 
ated by the central entrapped bubble, our images are 
always nearly uniform and dark, see second panel of Fig. 
Id. By carefully measuring the intensity in this region, 
we can constrain the height of any possible air film un- 
derneath the drop. 

For a liquid layer separated from a glass substrate by 
an air gap of height ft, the total electric field produced 
by multiple reflections from the two interfaces is: 
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where n g = 1.52 and n\ = 1.44 + 0.0032z are respec- 
tively the glass and liquid complex indicies of refraction. 
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The optical path length is given by S 
term e %5k in eqn. 2 accounts for the phase shift caused 
by the path-length difference 2hk after the k th reflection 



from the lamella. Only the first few terms in eqn. (2) are 
large enough to contribute significantly to the total elec- 
tric field. Thus, the finite coherence length of the LED 
does not influence this calculation. The total intensity 
is then = \El(K)\ 2 + i^, where 1^ is an unknown 

background intensity from stray light and incoherent re- 
flections. 

To obtain the air gap thickness, ft, we must eliminate 
the unknown quantities \Eq\ 2 and 1^ from our expression 
for II. To do this, we measure (i) the intensity under 
the lamella long after spreading has finished so that the 
liquid can be assumed to be in contact with the substrate, 
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and (ii) the intensity due to reflected light from only the 
substrate: 
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We can determine ft by computing: 

I L (h) - Ic 
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Fig. [3)2 shows I s from the time the cavity collapses 
into a bubble until just before sheet ejection. Fig. 3b-d 
shows I s measured just before the instant of sheet ejec- 
tion vs. i^o, v, P. In all cases, the data remain essentially 
constant and within error of zero. The error bars in Fig. 
[3] represent drop-to-drop fluctuations. The distribution 
of all measurements of I s has a mean 0.0004 with a stan- 
dard deviation 0.0054. This is comparable to the noise in 
the camera between adjacent frames when filming a still 
drop. 

These measurements of I s are consistent with the liq- 
uid being in direct contact with the substrate; this data 
places an upper bound of 3.8 nm (using eqn. [5]) on the 
thickness of any possible air layer beneath the spreading 
lamella for all of the parameter space sampled. An air 
layer of this thickness would be highly unstable and it is 
difficult to conceive that it could persist over 0.4 ms, i.e. 
until the moment of sheet ejection. All of the air trapped 
beneath the falling drop is enclosed into the small cen- 
tral bubble discussed above and does not influence the 
subsequent sheet ejection and splashing. 

However, once the thin sheet is ejected, it does move 
over a layer of air, which is easily visualized with our 
interference technique, see the bottom three panels in 
Fig. Id. The sheet is ejected at a very shallow angle, 
varying from 0.1° to 0.25°, and this angle is relatively 
insensitive to v and P but decreases with increasing uq. 

It is highly ant i- intuitive that the surrounding gas con- 
trols splashing in all viscosity regimes. Our interference 
technique allows quantitative measurements of the air be- 
neath a spreading drop as a function of position and time. 
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FIG. 3: Interference signal, I s , used to determine air-gap thickness, h, underneath the spreading lamella, (a) I s vs. time, from 
the time the cavity collapses until just before sheet ejection, (b-d) show I s measured one frame (~ 50 [is) before sheet ejection 
as a function of (b) impact velocity, (c) viscosity, and (d) pressure. These measurements place an upper bound on h of 3.8 nm. 
The data is consistent with the lamella making direct contact with the substrate as it expands. Error bars represent drop to 
drop fluctuations. 



In the low- viscosity splashing regime, corona formation 
occurs very near to the entrapped air bubble, both spa- 
tially and temporally. Techniques such as total internal 
reflection imaging have been used to explore low- viscosity 
impact dynamics [23 , confirming theoretical predictions 
[20] of the initial, transient air film. In higher viscosity 
fluids, we find these initial air-cavity dynamics are also 
in quantitative agreement those predictions. 

However, we find no trapped air beneath the spreading 
drop outside the small central bubble; there is no signif- 
icant air film beneath the drop at the time of thin-sheet 
ejection. This suggests that, rather than an underlying 
air layer, gas flow at the edge of the spreading drop is 
responsible for destabilizing the liquid. This conclusion 
is consistent with previous splash experiments in the low- 
viscosity regime [10] . In that case, the scaling with gas 
pressure and molecular weight suggests that the liquid 
front expanding into the surrounding air leads to liquid 
destabilization and splashing. The results reported here 
for more viscous fluids suggest a similar instability due 
to leading-edge gas flows. 
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